Direct observation of Dirac interface states at the Bi 2 Se3/Si(lll) interface 
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Topological insulators with mass-less helical Dirac fermions on their surface have attracted con- 
siderable attention in recent years. Several exotic effects and proposed applications rely on the 
presence of such states also at the interface between a topological insulator and a trivial insulator. 
In the present work, we have used low-energy angle-resolved photoelectron spectroscopy to demon- 
strate and characterize the interface states at the substrate interface in Bi2Se3 thin films grown on 
Si(lll). The results establish that Dirac fermions are indeed present also at the interface. It is also 
shown that the interface states are shifted in energy with respect to the surface states due to band 
bending in the film. The observation of topological interface states paves the way for further studies 
of interface related phenomena involving topological insulator thin films. 



The three-dimensional (3D) topological insulators 
(TIs) constitute a class of materials in which the inverted 
band structure of the insulating bulk is accompanied by 
the existence of conducting metallic states at the sur- 
face [1-3]. These surface states consist of massless he- 
lical Dirac fermions which are protected by time rever- 
sal symmetry and therefore robust against non-magnetic 
backscattering. Although discovered only recently, the 
TIs have been the subject of extensive research, not only 
due to their inherently interesting properties and po- 
tential impact on the development of room-temperature 
spintronic devices 0, [H| but also because they are pre- 
dicted to make the realization of a number of exotic phys- 
ical phenomena, such as Dyons jfjj|, Majorana fermions Q 
and axion dynamics [8|, feasible. The binary compound 
Bi2Se3, considered to be the model system for the 3D TI 
material class [9j, has received a lot of attention in this 
context. Photoemission studies of single crystal Bi2Se3 
have confirmed the nearly linear dispersion of the Dirac 
fermions [10] while the helicity has been investigated by 



recent spin-resolved measurements [U|, LL2| . 

Pioneering works on Bi2Se3 thin films have revealed a 
strong variation in the electronic band structure of the 
surface state as a function of film thickness [13| . At thick- 
nesses below 5 quintuple layers (QL), one QL is 9.5 A 
thick, an energy gap opens at the Dirac point as a result 
of hybridization between Dirac states originating from 
opposite sides of the film. The strength of the hybridiza- 
tion decreases rapidly with increasing film thickness and 
for 6 QL films it is negligible and the energy gap at the 
Dirac point closes. Experimental observations of this gap 
opening in ultra thin films provide indirect evidence for 
the existence of a Dirac state at the interface towards the 
substrate. In fact, theory predicts the existence of gapless 
topological Dirac states at interfaces between topologi- 
cally trivial materials and TIs [l4|. Consequently, such 
states should exist at the interface between a TI thin film 
and the substrate on which it is grown. However, so far 
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no direct observation of buried-interface states in non- 
hybridized TI films has been made. It is therefore, first 
of all, of fundamental interest to experimentally establish 
the existence of interface states in a Tl-semiconductor 
junction. Additionally, a qualitative description of how 
the electron configuration in thin films is influenced by 
the presence of the substrate, including a band alignment 
model for the junction, is desirable, something which thus 
far has only been discussed briefly 15 -I?)- Once under- 
stood, the latter can provide insights on how to tailor 
the electronic properties, e.g. the position of the chemi- 
cal potential, of a TI film by changing properties of the 
substrate. Such tailoring can prove useful not only in 
future realization of some of the mentioned exotic phe- 
nomena but also for practical device applications. In 
this Letter we address these issues and provide experi- 
mental evidence for the existence of Dirac fermions at 
the substrate interface in thin films and show how the 
band bending through the film can be explained by band 
alignment between the film and the substrate at the in- 
terface. 

In the present work we have studied the electronic band 
structure of Bi2Se3 thin films by means of laser-based 
Angle- Resolved Photoelectron Spectroscopy (ARPES). 
In order to enable direct observation of the buried- 
interface state a combination of low photon energy and 
cryogenic sample temperature was used in the experi- 
ments since this dramatically increases the escape depth 
of the photoexcited electrons [18]. The photoemission 
experiments were performed at the BALTAZAR laser- 
ARPES facility [Taj at a photon energy of 10.5 eV. Thin 
films of thicknesses 6 QL and 8 QL were grown in situ 
following the method described in Ref. [13] , in which thin 
films are deposited by co-evaporation of Bi and Se onto a 
Bi-terminated Si(lll)-(7 x 7) substrate. The specific film 
thicknesses were chosen so that hybridization between in- 
terface and surface Dirac states is avoided while, at the 
same time, the films are sufficiently thin for the interface 
to be directly probed by the ARPES experiment. 

Figure [T] presents photoemission data of the 8 QL and 
6 QL Bi2Se3 films. Both samples were grown on n-type 
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FIG. 1. (color online) (a)-(c) Photoemission intensity plots (raw data) of a 8 QL thick Bi 2 Se3 film on n-type Si-substrate taken 
at hv — 10.5 eV and T = 9 K. (a) Energy dispersion along the K-f-K direction. The fitted band structure described by Eq. 
[1] is plotted as solid lines on top of the experimental data, (b) Constant energy cut at Ef displaying the hexagonal Fermi 
surface, (c) Constant energy cut at i?b=0.25 eV. (d)-(f) Photoemission data of a 6 QL Bi2Se3 film on n-type Si-substrate 
taken at hv— 10.5 eV and T = 9 K. All energy distribution curves (EDCs) in the data set are normalized to unity in order 
to increase visibility of the interface state, (d) Measured energy dispersion along the K-f-K direction together with the fitted 
band structure of Eq. [2] Solid and dashed lines represent the surface and interface states, respectively, (e) Constant energy 
cut at Ef. Two hexagonal energy contours are observed, the inner one belonging to the surface state and the outer to the 
interface state, (f) Constant energy cut at E\>= 0.25 eV. The black areas in (e) and (f) are outside the detector, which explains 
why there are no counts in these areas. 



doped (0.004 Ocm) Si(lll) substrates and the data were 
taken at T = 9 K. The energy- momentum slice for the 8 
QL film, displayed in Fig. QJa), shows a nearly linearly 
dispersing Dirac surface state with the Dirac point lo- 
cated at a binding energy (E\>) of 0.46 eV. As revealed 
by the constant energy cut taken at the Fermi level (E-p), 
shown in Fig. DJb), the Fermi surface (FS) is hexagonally 
deformed. Closer to the Dirac point the constant energy 
surface becomes circular, as seen in Fig. Die). This de- 
viation from an ideal Dirac cone has been reported pre- 
viously [loj and is expected for highly n-doped Bi2Se3 
samples. By using the single band model [21] described 

by 



E±(k) = E -Dk 2 ± 




Bk 2 



(v F hk) 2 , (l) 



together with peak positions of momentum-distribution 
curves (MDCs) extracted from our data, we can create 
a fitted curve of the surface state. Here, E± represent 
the dispersion of the Dirac state above and below the 
Dirac point respectively, Eq is the binding energy of the 
Dirac point, A is the energy gap at the Dirac point (non- 
zero only if inter-surface coupling is present), is the 
Fermi velocity and k the in-plane electron momentum. 
B and D are model parameters and their definitions are 
found in Ref. |2l|. The resulting fitted curve is plot- 



ted on top of the experimental data in Fig. QJa) (solid 
line) and the fitted model parameters are listed in Ta- 
ble 1. The model described by Eq. (pQ) is derived from 
an effective Hamiltonian and fits well with the data for 
small k- values. However, it does not include the hexag- 



onal distortion of the Dirac cone. Therefore, only points 
in the region where no distortion of the cone is observed 
are included in the fit. Consequently, the resulting Fermi 
velocities (vp) from the fits, listed in Table HI are smaller 
than reported elsewhere OIlB- For comparison, actual 
values of vp are also given in Table HI 



Comparing the energy-momentum slices in Figs. QJa) 
and [l](d) for the 8 QL and 6 QL films respectively, re- 
veals the presence of similar Dirac-like states in both 
films with the Dirac points located at nearly the same 
binding energies. However, in the 6 QL data there is an 
additional, much larger, V-shaped feature reminiscent of 
a Dirac state shifted towards higher binding energy. The 
constant-energy slice displayed in Fig. [TJe) is taken at Ep 
and shows that both of these features are hexagonally de- 
formed similar to what we observe in the 8 QL film. For 
a constant-energy cut at higher binding energy the in- 
ner feature becomes circular, as seen in Fig. D](f), since 
we are now approaching the Dirac point of this state, 
while the outer constant-energy surface remains hexago- 
nal. The model fitted to the 8 QL data describes a single 
Dirac state only. Therefore, in order to determine the 
origin of the additional spectral feature seen in the 6 QL 
case we adopt the double band model from Ref. [17j . 
This model includes Dirac cones on opposite sides of the 
TI film and accounts for the structure inversion asymme- 
try (SI A) and the effective band bending induced by the 
presence of the substrate. The dispersion of the Dirac 
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FIG. 2. (color online) (a) Energy dispersion along the K-f-K 
direction for a 6 QL Bi2Se3 grown on a n-type substrate, (b) 
Corresponding energy- momentum slice for a 6 QL film on p- 
type substrate. Due to missing data for k y values below -0.25 
A -1 the intensity in this region has been manually set to zero. 



states is then given by 




E a± (k) = E - Dk 2 ±\l(j- Bk 2 ) + (\V\ + av F mm 

The term \ V\ represents the band bending across the film 
and together with a = ±1 it produces a set of two Dirac 
cones with opposite helicities having their Dirac points 
separated in energy by 2\V\. Eo is now the center energy 
of the two Dirac points. The remaining parameters have 
the same meaning as in Eq. (pp). 

Fitting Eq. ((2]) to our data reproduces the observed 
band structure well, as seen in Fig. QJd). Also in this 
case, only data points close to the Dirac points are used 
in the fit. The large V-shaped state observed in the sim- 
ulated band structure (dashed line in Fig. DJd)) corre- 
sponds to the Dirac cone at the interface between the TI 
and the substrate, which suggest that what we observe 
is indeed the interface state of the TI film. The fact that 
no such state is observed in the 8 QL ARPES data in- 
dicates that the thickness of this film is larger than the 
electron escape depth at this particular excitation energy 
and temperature. 

In order to verify the interface nature of the observed 
feature in the 6 QL film we have performed photon en- 
ergy dependent ARPES measurements at the 13 beamline 
at MAXlab [22| . A 6 QL film was measured at a photon 
energy of 10 eV displaying both the surface and the inter- 
face states. However, when increasing the photon energy 
to 22 eV only the surface state remained visible. The re- 
duced intensity of the interface state in the photoelectron 
spectrum for higher excitation energy is consistent with 
the increased surface sensitivity of the ARPES technique 
in that particular photon energy range. 

All together, the observed thickness and photon energy 
dependence of the states in the 6 QL film, combined with 
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TABLE I. Fitted parameters for the 6 QL films using Eq. © 
and for the 8 QL film using Eq. (pQ). Actual values for vf 
along the f-K direction, based on MDC peaks close to Ef, 
are also listed. 



the agreement between the theoretical model and experi- 
mental data, support the notion that we directly observe 
the Dirac state located at the substrate interface. The 
interface state exhibits a similar Dirac-like dispersion as 
the surface state and is shifted towards higher binding 
energy as a result of band bending in the film. 

In order to gain a qualitative understanding of the sub- 
strate effect, and hence the band bending in the film, we 
have studied 6 QL films grown on n-type (0.004 Ocm) 
and p-type (0.0009 Qcm) Si(lll) substrates. Comparison 
reveals, interestingly, very little differences between the 
observed band structure of the films, see Fig. [2j although 
the bulk Fermi level (Ev) position in the two substrates 
are very different. In both cases the energy separation 
between the two Dirac points (A .Ed) is roughly 0.3 eV. 

The striking similarities between the 6 QL films on n- 
type and p-type substrates is a consequence of the Fermi 
level of the Si(lll)-(7 x 7) surface being pinned approx- 
imately 0.7 eV above the valence band maximum (Ey) 
due to a high density of states located in the band gap at 
the surface [23l425j. This pinning of E F is independent 
of bulk doping and thus equivalent for both n-type and 
p-type substrates. Terminating the surface with 1 ML 
bismuth has little effect on the pinning level, reducing 
E F — E v from 0.7 eV to approximately 0.65 eV [26[. 

From the measured value of AEd we get the magni- 
tude of the band bending across the film. The binding 
energy positions of the Dirac points, together with the 
known band gap of Bi 2 Se 3 (0.35 eV [jUSzl), 

gives the 

Fermi level relative the conduction band minimum. As- 
suming an unchanged position of Ev at the substrate 
surface, alignment of the Fermi levels in the film and 
substrate results in conduction band offsets (AEc) of 
1.09 eV and 1.02 eV for the films on n-type and p-type 
substrate respectively. Here we have also assumed that 
the Dirac point is located in the center of the Bi2Se3 
band gap. Figures [31(a) and [31(b) display the resulting 
band alignment for the Bi2Se3/Si heterojunction. The 
upwards (downwards) band bending seen for the n-type 
(p-type) substrate is a result of the discussed pinning 
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FIG. 3. (color online) (a) Band diagram for the Bi 2 Se3/Bi/n- 
Si system showing the conduction (Ec) and valence (Ey) 
bands in the thin film and substrate. The conduction band 
offset AEc at the interface is 1.09 eV. The dashed horizon- 
tal line refers to the Fermi level, Ef> Vertical lines represent 
the vacuum/TI and Tl/substrate interfaces respectively. A 
schematic drawing of the observed band structure for the sur- 
face and interface states with corresponding values for En 
is enclosed by the dashed rectangle, (b) Band diagram for 
the Bi 2 Se 3 /Bi/p-Si system for which A£ c =1.02 eV. (c) Pre- 
dicted band diagram for a p-doped Bi2Se3 thin film on a H- 
terminated Si(lll) substrate. 



of E-p at the surface. Therefore, the energy difference 
E-p — E v at the substrate side of the interface is the same 
in the two cases (0.65 eV). For simplicity, the band bend- 
ing through the thin film is assumed to be linear. If we 
instead use the electron affinity rule [28] to determine the 
conduction band offset we obtain AEq = 1.14 eV. The 
electron affinities x[Si(lll)-(7 x 7)]= 4.16 eV [29] and 
x[Bi2Se3]~5.3 eV have been used. The latter is an esti- 
mate based on the value of the work function of Bi2Se3 
which is determined from the photoemission experiment. 
The similar values for AEc obtained using the two ap- 
proaches indicate that the sketched band alignment is 
qualitatively correct and that the pinned Fermi level po- 
sition of the substrate is little influenced by the presence 
of the film. Also, the magnitude of AEc at the Bi2Se3-Si 
interface is comparable to the band gap of silicon. 

As seen from our results, using Si(lll) substrates with 
very different doping levels does not influence the band 
bending through the Bi2Se3 thin films, as initially ex- 
pected, due to the pinning of the Fermi level at the in- 



terface. This seems to limit the possibilities of exploit- 
ing the bulk position of Ep in the substrate to modify 
the electronic configuration of TI thin films. However, 
this problem might be overcome by instead using a H- 
terminated substrate. At the H-Si(lll) surface the pin- 
ning of Ep is removed 30|, 3l[ and the unpinned position 
is determined by the Fermi level in the bulk. Tuning the 
position of Ep in the substrate by bulk doping can then 
be used to manipulate the band bending across the film. 

Since the conduction band offset AEq at the interface 
is fixed, tuning of the chemical potential in the film by 
Ca or Mg-doping [32|, [33| is expected to shift the surface 
Dirac point towards lower binding energy. We propose 
the idea of using a heavily p-doped H-Si(lll) substrate 
combined with a p-doped Bi2Se3 film to obtain the band 
alignment presented in Fig. [3fc). If the film is suffi- 
ciently doped, the chemical potential at the surface can 
be placed below the Dirac point which would result in a 
hole-like surface state and electron-like interface state as 
shown in Fig. [3^c). This system could then host bound 
electron-hole pairs, or excitons, and possibly permit the 
observation of a topological exciton condensate (TEC) as 
suggested in Ref. [34|. 

The work presented in this article establishes exper- 
imentally the existence of interface states at the inter- 
face between a trivial insulator and a topological insula- 
tor. Similar to the surface state on Bi2Se3 the interface 
state exhibits a Dirac-like linear dispersion. The band 
bending through the film, generated by the substrate, 
shifts the Dirac point of the interface state in energy rel- 
ative the surface state. We predict that H-termination 
of the Si(lll) substrate will allow the band bending 
across the film to be manipulated by changing the dop- 
ing level and carrier type of the substrate. Our work also 
shows that the combination of low-temperature and low- 
photon-energy ARPES permits direct studies of buried 
interfaces in these materials, something which opens up 
for further studies on, for example, interfaces in p-n TI 
junctions and Tl-heterojunctions. 

This work was made possible through support from the 
Knut and Alice Wallenberg Foundation and the Swedish 
Research Council. 



[1] L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 
106803 (2007). 

[2] J. E. Moore and L. Balents, Phys. Rev. B 75, 121306(R) 
(2007). 

[3] L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 
[4] O. V. Yazyev, J. E. Moore, and S. G. Louie, Phys. Rev. 

Lett. 105, 266806 (2010). 
[5] H. Steinberg, J.-B. Laloe, V. Fatemi, J. S. Moodera, and 

P. Jarillo-Herrero, Phys. Rev. B 84, 233101 (2011). 
[6] X.-L. Qi, R. Li, J. Zang, and S.-C. Zhang, Science 323, 

1184 (2009). 



5 



[7] L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 

(2008) . 

[8] R. Li, J. Wang, X.-L. Qi, and S.-C. Zhang, Nature Phys. 
6, 284 (2010). 

[9] H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and 

S.-C. Zhang, Nature Phys. 5, 438 (2009). 
[10] Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, 

A. Bansil, D. Grauer, Y. S. Hor, R. J. Cava, and M. Z. 

Hasan, Nature Phys. 5, 398 (2009). 
[11] Z.-H. Pan, E. Vescovo, A. V. Fedorov, D. Gardner, Y. S. 

Lee, S. Chu, G. D. Gu, and T. Valla, Phys. Rev. Lett. 

106, 257004 (2011). 
[12] C. Jozwiak, Y. L. Chen, A. V. Fedorov, J. G. Analytis, 

C. R. Rotundu, A. K. Schmid, J. D. Denlinger, Y.-D. 

Chuang, D.-H. Lee, I. R. Fisher, R. J. Birgeneau, Z.-X. 

Shen, Z. Hussain, and A. Lanzara, Phys. Rev. B 84, 

165113 (2011). 

[13] G. Zhang, H. Qin, J. Teng, J. Guo, Q. Guo, X. Dai, 
Z. Fang, and K. Wu, Appl. Phys. Lett. 95, 053114 

(2009) . 

[14] J.-H. Song, H. Jin, and A. J. Freeman, Phys. Rev. Lett. 

105, 096403 (2010). 
[15] Y. Zhang, K. He, C.-Z. Chang, C.-L. Song, L.-L. Wang, 

X. Chen, J.-F. Jia, Z. Fang, X. Dai, W.-Y. Shan, S.-Q. 

Shen, Q. Niu, X.-L. Qi, S.-C. Zhang, X.-C. Ma, and 

Q.-K. Xue, Nature Phys. 6, 584 (2010). 
[16] Y. Sakamoto, T. Hirahara, H. Miyazaki, S. I. Kimura, 

and S. Hasegawa, Phys. Rev. B 81, 165432 (2010). 
[17] W.-Y. Shan, H.-Z. Lu, and S.-Q. Shen, New J. Phys. 12, 

043048 (2010). 

[18] M. P. Seah and W. A. Dench, Surf. Interface Anal. 1, 2 
(1979). 

[19] M. H. Berntsen, O. Gotberg, and O. Tjernberg, Rev. 
Sci. Instrum. 82, 095113 (2011). 

[20] K. Kuroda, M. Arita, K. Miyamoto, M. Ye, J. Jiang, 
A. Kimura, E. E. Krasovskii, E. V. Chulkov, H. Iwasawa, 
T. Okuda, K. Shimada, Y. Ueda, H. Namatame, and 
M. Taniguchi, Phys. Rev. Lett. 105, 076802 (2010). 



[21] H.-Z. Lu, W.-Y. Shan, W. Yao, Q. Niu, and S.-Q. Shen, 
Phys. Rev. B 81, 115407 (2010). 

[22] M. H. Berntsen, P. Palmgren, M. Leandersson, A. Hahlin, 
J. Ahlund, B. Wannberg, M. Mansson, and O. Tjern- 
berg, Rev. Sci. Instrum. 81, 035104 (2010). 

[23] F. J. Himpsel, D. E. Eastman, P. Heimann, B. Reihl, 

C. W. White, and D. M. Zehner, Phys. Rev. B 24, 1120 
(1981). 

[24] F. J. Himpsel, G. Hollinger, and R. A. Pollak, Phys. 

Rev. B 28, 7014 (1983). 
[25] J. M. Nicholls and B. Reihl, Phys. Rev. B 36, 8071 

(1987). 

[26] K. Hricovini, G. Le Lay, A. Kahn, A. Taleb-Ibrahimi, 

J. E. Bonnet, L. Lassabatere, and M. Dumas, Surf. Sci. 

252, 424 (1991). 
[27] J. Black, E. M. Conwell, L. Seigle, and C. W. Spencer, 

J. Phys. Chem. Solids 2, 240 (1957). 
[28] R. L. Anderson, Solid-State Electron. 5, 341 (1962). 
[29] G. Hollinger and F. J. Himpsel, J. Vac. Sci. Technol. A 

1, 640 (1983). 

[30] C. J. Karlsson, E. Landemark, L. S. O. Johansson, U. O. 
Karlsson, and R. I. G. Uhrberg, Phys. Rev. B 41, 1521 
(1990). 

[31] R. Hunger, C. Pettenkofer, and R. Scheer, J. Appl. Phys. 
91, 6560 (2002). 

[32] D. Hsieh, Y. Xia, D. Qian, L. Wray, J. H. Dil, F. Meier, 
J. Osterwalder, L. Patthey, J. G. Checkelsky, N. P. Ong, 
A. V. Fedorov, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, 
R. J. Cava, and M. Z. Hasan, Nature (London) 460, 
1101 (2009). 

[33] Y. L. Chen, J.-H. Chu, J. G. Analytis, Z. K. Liu, 
K. Igarashi, H.-H. Kuo, X. L. Qi, S. K. Mo, R. G. Moore, 

D. H. Lu, M. Hashimoto, T. Sasagawa, S. C. Zhang, I. R. 
Fisher, Z. Hussain, and Z. X. Shen, Science 329, 659 
(2010). 

[34] B. Seradjeh, J. E. Moore, and M. Franz, Phys. Rev. Lett. 
103, 066402 (2009). 



